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^ \ Abstract 

-Y-\ . To describe the dynamics of a single peak of the Rosensweig instabihty a model is 

proposed which approximates the peak by a half-ellipsoid atop a layer of magnetic 
^^ ■ fluid. The resulting nonlinear equation for the height of the peak leads to the correct 

*yN . subcritical character of the bifurcation for static induction. For a time-dependent in- 

<^ \ duction the effects of inertia and damping are incorporated. The results of the model 

show qualitative agreement with the experimental findings, as in the appearance of 
period doubling, trebling, and higher multiples of the driving period. Furthermore 
^S \ a quantitative agreement is also found for the parameter ranges of frequency and 

induction in which these phenomena occur. 
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1 Introduction 



Magnetic fluids (MF) are stable colloidal suspensions of ferromagnetic nanopar- 
ticles (typically magnetite or cobalt) dispersed in a carrier liquid (typically oil 
or water). The nanoparticles are coated with a layer of chemically adsorbed 
surfactants to avoid agglomeration. The behaviour of MF is characterized by 
the complex interaction of their hydrodynamic and magnetic properties with 
external forces. Magnetic fluids have a wide range of applications [1] and show 
many fascinating effects [2], as the labyrinthine instability or the Rosensweig 
instability. The latter instability occurs when a layer of MF with a free surface 
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is subjected to a uniform and vertically oriented magnetic field. Above a cer- 
tain threshold of the magnetic field that surface becomes unstable, giving rise 
to a hexagonal pattern of peaks [3] . Superimposing the static magnetic field 
with oscillating external forces leads to nonlinear surface oscillations. Exper- 
imentally, either vertical vibrations [4-6] or magnetic fields [7-11] have been 
investigated as alternating external forces. 

For free surface phenomena the fluid motion strongly depends on the shape 
of the surface and vice versa. Additionally for MF, the shape of the surface 
is determined by the magnetic field configuration which contributes via the 
Kelvin force to the Navier-Stokes equation the solution of which gives the 
flow field. Thus the dynamics of MF is inherently governed by the nonlinear 
interaction between the fiow field, the surface shape, and the magnetic field 
configuration. 

For that reason one attempts to study simple systems of MF which never- 
theless show the essential features. The nonlinear dynamics of a single peak 
of magnetic fiuid, i.e., the dynamics of a 0-dimensional system in a vertically 
oscillating magnetic field was studied exemplarily in [10]. By varying the am- 
plitude and the frequency of the alternating field and the strength of the static 
field, the peak response can be harmonic, subharmonic (twice the driving pe- 
riod) or higher multiples of the driving period. For suitable choices of the 
parameters, non-periodic chaotic peak oscillations were observed. 

Taking the above described circumstances into account for a theoretical ap- 
proach, a sound model should be analytically tractable as well as capable of 
showing all essential features. Beyond these primary demands, the model may 
also predict new phenomena of peak oscillations. The aspiration to confirm 
such new phenomena experimentally motivates a simple and robust model to 
guide the design of the experimental setup. 

Such a model is proposed for the dynamics of a single peak of MF. It is based 
on the approximation of the peak by a half-ellipsoid with the same height 
and radius as the peak. The resulting equation giving the dependence of the 
height of the peak on the applied induction is derived in the following section. 
The character of the bifurcation is analysed in Sec. Ill for the case of a static 
induction. In Sec. IV the dynamics of the peak is studied and the results are 
compared with the experimental behaviour for different frequencies of a time- 
dependent induction. The final section summarizes the results and outlines 
two aspects for further experiments. 



2 Model 



The complex and nonlinear interactions in MF with a free surface formed by 
a peak (see Fig. 1) present a formidable problem, since the form of the peak 
is not known analytically. The aim of our model is an analytical equation for 
the height of the peak at its centre r = {x, y) = 0. The shape of the peak, 
particular the form at the tip of the peak is beyond the potential of this model. 
The equation will thus neglect the influence of the surface regions away from 
the peak tip and of the boundaries. 

A layer of an incompressible, nonconducting, and inviscid MF of half-infinite 
thickness between z = and z -^ — oo is considered with a free surface 
described hj z = ({x, y, t). It is assumed that the magnetization M of the MF 
depends linearly on the applied magnetic field H, M = ^H, where x is the 
susceptibility of the MF. The system is governed by the equation of continuity, 
divv = 0, and the Euler equation for MF in the presence of gravity 



p [dtv + (vgrad) v] = — gradp + /iq MgradiJ 
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where the magnetostriction is neglected and the co-linearity of the magneti- 
zation and the field is exploited for the magnetic force term. In Eq. (1) the 
velocity field is denoted by v, the density of the MF by p, the pressure by p, 
the permeability of free space by po and the acceleration due to gravity by g. 
M, H, and B are the absolute values of the magnetization, the magnetic field 
and the induction B in the fluid. In the static case, v = 0, the integral of the 
equation of motion (1) may be calculated to give 

H 

p=-pgz + poMdH' + const. (2) 





The remaining boundary condition in the static case, the continuity of the 
normal stress across the free surface, gives 

p = aK-^^{Mnf e.t z = ( , (3) 



where the pressure in the non-magnetic medium above the MF was set to 
zero. The surface tension between the magnetic and non-magnetic medium is 
denoted by a, the curvature of the surface hj K = div n, and the unit vector 
normal to the surface by n. By inserting Eq. (2) at z = ( into Eq. (3), the 



balance of pressure results in an equation for the surface ( 
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with M„ = Mn. After the peak is formed, the equilibrium is characterized 
by the equality of the pressure along the surface. Motivated by our aim of an 
analytically tractable model, we choose the two reference points r = 0, the 
centre of the peak, and |r| ^ 0, the fiat interface far away from the peak, 
where the pressure equahty is evaluated. The magnetization is related to the 
induction by 
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Applying Eq. (4) at (r = 0, C = h) and (|r| > 0, C = 0) leads to 
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where -Bext is the external applied induction. The remaining two unknown 
quantities, the curvature K{h) and the induction B{h) at the tip of the peak, 
are determined by an approximation. We model the peak, which is assumed to 
be rotationally symmetric, by a half-ellipsoid with the same height and radius 
as the peak (see Fig. 1). Thus one can make use of the analytical results for 
a rotational ellipsoid with the vertical (horizontal) semiprincipal axis h (R) 
with the curvature given by 
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and the induction [12] 



P={ 



l + e' 



l-e' 



^ e- 



B 



B. 



ell 



(e — arctanel 



(artanh 



e -e 



X+1 

i + x/3 



Bp^t with 



./im^^i 



v/i - Why 



ioT R> h 



(9) 



ioT R < h. 



It is emphasized that an apphed induction Bg^t results in a uniform induction 
Bell within the elhpsoid. The demagnetization factor /9 is a purely geometrical 
quantity because it relates the dimensions of the major and minor semiprinci- 
pal axes by means of the eccentricity e. Whereas (7) can be substituted directly 
into Eq. (6), the result (8) has to be modified to the case of a half-ellipsoid 
atop the layer of MF. The proposed modification is 
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where a parameter A is introduced, which mimics the influence of the mag- 
netic field of the layer on the field at the tip of the peak. The form of (10) 
ensures that in the limits of a magnetically impermeable material {x = 0); 
of a 'magnetic conductor' (x -^ oo), of a very oblate ellipsoid (/5 — 1), and 
of a very prolate ellipsoid {j3 ^ 0) the results are the same as in Eq. (8). As 
long as the height of the half-ellipsoid is large compared to its diameter, the 
influence of the magnetic layer on the magnetic field at the tip of the peak is 
small. This is obviously not the case if the half-ellipsoid becomes disk-shaped, 
i.e. h < R. For this case (10) is expanded up to the first order in h/R, 
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and compared to the analytical result for B at the crest of a sinusoidal surface 
wave (SW) (pp. 178 in [2]) with the wave length 4R (Fig. 2) 
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The condition that both values of B should coincide, leads to an equation for 
the parameter A 
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The determination of A adjusts the radius, since the critical wave number for 
surface waves is the capillary wave number, kc = {2n/\c) = J pg/a. Therefore 
the radius of the half-ellipsoid is fixed to R = (Ac/4) = 7r/(2 k^). By inserting 
(7) into (6) and introducing dimensionless quantities for all lengths and the 
induction 
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we obtain a nonlinear equation for the dependence of the peak height h on 
the apphed induction Bg^t (the bars are omitted for the rest of the paper) 



r2 

^ext 




h 



1 



^ 0. (15) 



(x + 2) 



The nonhnear behaviour enters into the equation through B{h)/ B^xt which is 
determined by (9, 10, 13). Eq. (15) presents the fundamental equation of the 
model in which the height of the peak depends on the properties of the applied 
induction only. The quality of the approximation is tested in the static case 
for which (15) was derived. It forms the starting point for the description of 
the peak dynamics, where the effects of inertia and damping have to be taken 
into account. 



3 Static Peak 



For a layer of MF with a free surface subjected to a vertical magnetic field 
there are three different energies which contribute to the total energy Etot- 
The potential energy and the surface energy increase Etot with increasing h, 
whereas the magnetic field energy decreases Etot with increasing h. The plane 
surface corresponds to a minimum of the total energy at /i = 0. If the surface is 
perturbed, the magnetic flux is concentrated in the peaks of the disturbances. 
The resulting force tends to increase the modulations, while surface tension 
and gravitational forces tend to decrease the disturbances. When the increas- 
ing field passes a certain strength He-, the destabilizing force will win over 
the stabilizing ones. The resulting peaks are energetically favourable because 
for H > He the total energy has now a second minimum at /i > which is 
deeper than the first one at h = 0. The transition from the first to the second 
minimum corresponds to the sudden jump from h = to h > 0. li the peaks 
are established, a decreasing field results in smaller heights of the peak up to a 
second critical field Hs, the saddle- node field, where the peaks suddenly break 
down. With respect to the total energy this means a transition back to the 
first minimum at /i = because it is now energetically more favourable. 

Such a dependence of the height of the peak on the variation of the magnetic 
field is typified as a hysteresis. The difference between the two critical fields 
defines the width of the hysteresis. For a MF with x = 1-15 the width was 
measured to 6% of the critical field and the critical height at He is given by 
2.1 mm [10]. The corresponding critical inductions are Be = Ho[Hc + M{Hc)] 
andBs = fio[Hs + M{H,)]. 

For a static induction, B^xt = Bq, the solution of Eq. (15) is determined for two 



susceptibilities, x = 1-15 and x = 2.5. The former value is given in [10] for a 
mixture of EMG 901 and EMG 909 (both Ferrofluidics Corporation) in a ratio 
of 7 to 3, whereas the latter value was measured in a recent experiment for 
the same mixture [13]. For both susceptibilities a distinct hysteresis appears, 
whose width increases with increasing susceptibility. Correspondingly, in the 
limit X ^ the hysteresis disappears (Fig. 3a). For x = 1-15 the width is 5% of 
the critical induction B^ and the critical height of the peak is he — 2.0/ kc — 2.9 
mm. The material parameters p = 1377 kgm~^ and a = 2.86 ■ 10~^ kgs~^ as 
given in [10] were used. For the other chosen susceptibility, x = 2.5, the width 
is 13% of the critical induction B^ and the critical height is he — 6.9/ kc — 10.0 
mm. 

Figure 3b shows the dependence of the width of the hysteresis on the suscep- 
tibility of the magnetic fluid. Whereas for small susceptibilities a fair increase 
of the width can be detected, a tendency towards a saturation in the growth 
can be seen for larger susceptibilities. No systematic measurements of the 
width of the hysteresis have yet been undertaken. Therefore any experimental 
test which would determine the range of validity of the model is pending on 
subsequent measurements. 

Despite the simplicity of the proposed approximation, the model describes the 
generic static behaviour of the height of the peak very well, i.e. the appearance 
of a hysteresis for increasing and decreasing induction at nonzero susceptibili- 
ties as it is observed in experiments [10,14,15]. Note in this connection that for 
A = 0, i.e. when neglecting the difference in the magnetic field of an ellipsoid 
and a half-ellipsoid, no hysteresis is found. Note also that in a one-dimensional 
system one finds a supercritical bifurcation for x < 2.53 [16,17] whereas our 
simplified two-dimensional model yields a subcritical bifurcation for all val- 
ues of X ill accordance with experiment. Beyond the qualitative agreement, 
the quantitative values for the width of the hysteresis and the critical height 
are in satisfying agreement with the measurements in [10] for x = 1-15. This 
agreement is achieved without any fit-parameter since the fixed value of the 
parameter A applies for any MF. 

The fact that the critical induction is not equal to one (cf. 15) is a consequence 
of the evaluation of the introduced parameter A, which determines the radius. 
The imposed value of the radius ensures the equality of the magnetic induction 
at the top of the oblate ellipsoid and the crest of the surface wave. But the 
curvature is different: h/R"^ at the top of the ellipsoid is smaller than the value 
7r^/i/(4i?^) it takes at the crest of the surface wave. Thus the expansion of (15) 
for small h with R = 7r/2 
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leads to a critical induction smaller than 1, Be = Jl/2 + 2/n'^ ~ 0.84. The 
half-ellipsoid approximation (10) with A = —1/2 was quantitatively compared 
to a numerically exact determination of the magnetic field of a rotational half- 
ellipsoid atop of a horizontal layer by solving the Laplace equation for the 
magnetic potential (Fig. 4). For 0.5 < h/R < 6.5 the magnetic induction at 
the tip of the peak is approximated with an accuracy of 1.7%. The comparison 
shows that the modification of the magnetic field at the tip of the peak through 
the magnetic field of the layer is rather weak even for small heights. This 
supports our assumption that the field at the tip of the peak is the essential 
feature to describe its behaviour. 

Therefore Eq. (10) describes B directly at the height of the peak fairly accu- 
rately. Furthermore, equation (15) leads to the correct character of the bifur- 
cation, i.e. a subcritical instability, and gives the right width of the hysteresis 
compared with the experimental results [10]. With this level of confirmation, 
the dynamics of a single peak of MF is studied in the following section. 



4 Oscillating Peak 

4-1 Inertia and Damping 



The induction is chosen to be a superposition of a static part, Bq, and a time- 
dependent part, ABcos{ujt). The amplitude of the oscillating part is denoted 
by AB and the frequency by cu = 27r/ = 2tt/T. In correspondence with the 
experiments [10], the response-period of the peak is studied in dependence 
on the three parameters, the strength of the static part, the amplitude of the 
alternating part, and the driving frequency. If the last two parameters are kept 
constant, one distinguishes between three different regimes for the behaviour 
of h{t) with increasing Bq. For small Bq the surface remains fiat, i.e. h{t) = 0. 
Beyond a first, lower threshold h{t) oscillates between zero and a maximum 
hmax whereas beyond a second, higher threshold it alternates between two 
positive extrema, < hmin < hmax (see Fig. 5). The behaviour in the second 
regime will be the focus of our study since it was analysed experimentally in 
detail in [10]. 

In order to formulate a differential equation for the peak dynamics, the effects 
of inertia and damping have to be incorporated into Eq. (15). Since each term 
in (15) stems from the equation of pressure balance, the inertial term may be 
written as 

force md?h p\h\AdPh \<\^'^h 



The sign of proportionality indicates that in the frame of our model the mass 
and the area of the peak cannot be precisely determined. For these quantities 
the knowledge of the complete surface and the flow field are necessary. For 
this reason we choose the simple relation of a linear dependence of the mass 
of the peak on its height. 

The implementation of the damping is difficult. In the experiment one observes 
that the peak periodically arises up to a maximal height and collapses to zero 
height. This behaviour leads to the assumption that the system is endowed 
with a dissipation mechanism which acts particularly strongly when the col- 
lapsing peak approaches z = 0. Since such a mechanism cannot be derived in 
the frame of the present model, the idea of an impact oscillator [19,20] is used. 
The impact oscillator is an externally excited oscillator, where the oscillating 
mass impacts on a fixed boundary. From this boundary the mass is refiected 
with a velocity 
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where r is the coefficient of restitution and to is the time of the impact, /i(to) = 
0. Consequently, there are oscillations between < h < oo only. For a weakly 
damped impact oscillator it is known [19] that infinite series of transitions 
from period 1 to period A^ (A^ = 3, 4, . . .) can appear. This phenomenological 
resemblance to the observations in [10] also motivates the use of the idea of 
an impact oscillator. It is emphasized that the chosen special form of damping 
applies only to the second regime, where h{t) oscillates between zero and a 
maximum hmax- 

In our model an impact with z = occurs whenever h{t) reaches zero. The 
height and the velocity after the impact are fixed and independent of the 
behaviour before the impact. We choose 

dh 
h = 10"^ and — = at t = V , (19) 

dt 



which corresponds to a nearly complete dissipation of the energy at every 
impact. A similar choice was made for a model proposed in [10]. The choice 
of fixed values is obviously an oversimplification because it does not make any 
difference whether a large peak with a high velocity rushes towards z = or 
whether a small peak slowly approaches z = 0. 

The resulting differential equation for this cut-off mechanism in dimensionless 



quantities is (the time is scaled by g^^^p^/'^a ^/^) 
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with Bext = -Bo + ABcos{ujt). Eq. (20) is solved by means of the forth-order 
Runge-Kutta integration method with a standard time step of T/200000. The 
other standard parameters for the integration are h{0) = 10~^ and dth{0) = 
as initial conditions and a total time of 200 T over which the solution is 
calculated. The first 100 periods are considered as transient time for the system 
to relax to a response-behaviour independent of the initial conditions. The last 
100 periods are analysed with respect to a periodic behaviour of h{t) by means 
of a Poincare section. For our one-dimensional dynamics a Poincare section 
means to compare /i at a certain time, say tm = i^T, with h at times, which 
are N periods (A^ = 1, 2, . . .) later with respect to tm'- 



h{tm) = h{tm + T) 
h(tm) = h(tm + 2T) 

h{t^) = h{tm + NT). (21) 

Those equations which are fulfilled give the period A^ (and any higher multiples 
of A^) of the peak response. The chosen 100 periods of analysis ensure a good 
reliability of the estimated periods up to 30. 



4-2 Results and Discussion 



The results of the Poincare sections are plotted as period diagrams in the 
Bq-AB plane at a fixed frequency / and for two different susceptibilities 
(see Figs. 6, 7, 9). The constant part, Bq, is sampled in steps of 0.01. The 
amplitude of the alternating part, AB, is increased in steps of 0.025 with an 
initial value of 0.05. The different periods in the interesting second regime, 
< h(t) < hmax, are coded by colours. The periods 1 to 10 are encoded 
by a chart of distinctive colours starting with green, red, blue, and ending 
with orange. The higher periods from 11 to 30 are encoded by a continuous 
colour chart. Periods above 30 or a non-periodic behaviour of h{t) are noted by 
grey. This selection of colours is guided by the choice of colours in [10]. White 
areas inside and right of the coloured horizontal stripes indicate regions, where 
h{t) oscillates between two positive extrema. White areas left of the coloured 
stripes denote the regime h{t) = 0. 
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The period diagram for a low frequency of / = 0.03 (^ 2.5 Hz) is displayed in 
Fig. 6. In accordance with the experimental results for 2.5 Hz (see Fig. 5a in 
[10]) the response of the peak is harmonic almost everywhere in the Bq-AB 
plane. Responses with higher periods are detected only at the edge towards 
the third regime. The area of harmonic response is cone-like shaped, where 
the limit to the left is given by AB = Be — Bo for Bq < Be (solid line) and the 
limit to the right is given by AB = Bo — Be for Bq > Be (dashed line). These 
strict limits apply particularly to the MF with x = 2.5 (Fig. 6b), whereas the 
right limit is more frayed for the MF with x = 1.15 (Fig. 6a). The feature of a 
cone-like shape is also found in the experiment, but with a slight asymmetry 
at very small amplitudes of the alternating field. An asymmetry could not be 
found with our simple model. Another difference is that our right limit is too 
large compared to the experimental data. 

The appearance of only harmonic responses is caused by the low frequency. 
The corresponding characteristic time of the excitation, T, is large enough 
for the peak to follow the slow modulations of the external field. Therefore 
the peak oscillates with the same frequency as the external excitation. By 
considering the quasi-static limit / — > (see Fig. 3), the boundaries of the 
second regime can be understood as follows: As long as Bq is smaller than 
Be — AB, the resulting external induction varies over a range where at its 
lower bound /ii = is stable. At its upper bound either /i2 = is stable if 
Bq + AB < Bs, or h2 = and /i2 > are stable if Bq + AB > Bg. Because 
of the greater attraction of zero height in the bistable area due to the strong 
damping at /i = 0, the dynamics of h(t) is bounded by zero in both cases. 
If Bq is larger than Be + AB, the resulting external induction varies over 
a range where at its lower bound /ii > is stable and at its upper bound 
/i2 > hi is stable. Thus the dynamics of h{t) is bounded between hi and /i2. 
Consequently, for low frequencies the peak alternates in the second regime as 
long as Be — AB < Bq < Be + AB. The fact that h{t) remains at zero even 
when for a certain time a nonzero height is stable (but not attractive enough 
to win over h = 0) was observed in the experiment, too. In our dynamics the 
zero height is always more attractive than the nonzero height. This is not the 
case in the experiment, which explains the observed lower limit to the right. 

Fig. 7 shows the results for a medium frequency of / = 0.1 (~ 8.2 Hz). For 
the MF with the low susceptibility of x = 1-15, the second regime splits into 
two disjoint parts. For smaller values of Bo we find only harmonic responses, 
whereas for higher values of Bo we observe the period N = 2. The second 
regime is separated from the first regime by AB = Be — Bq for Bq < Be 
(solid line). The limit to the right is given by AB = Bo — Be for Bq > Be 
(dashed line) only for amplitudes AB above 0.15 (Fig. 7a). For the MF with 
the high susceptibility of x = 2.5, the second regime forms a compact region, 
which is separated from the first regime by AB = Be — Bo for Bq < Be 
(solid line). In contrast to the low frequency behaviour, the whole structure 
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of periods shows a specific composition. For a fixed amplitude AB the peak 
starts to oscillate harmonically. For AB > 0.25 and increasing Bq the period-1 
state is replaced by the period-2 state which lasts up to the right limit of the 
second regime. This clear two-state picture changes if AB is decreased. For 
AB < 0.25 a tongue of high periodic (A^ > 2) and non-period oscillations 
appears (Fig. 7b). For 0.175 < AB < 0.25 the tongue is embedded in the 
period-2 state. For AB < 0.175 the tongue follows directly the harmonic 
oscillations. In this tongue we find odd-number periods of 3, 5, and 9 and 
even-number periods of 4, 8, 14, 16, and 18 (see Fig. 8). 

The structure of periods in Fig. 7b displays generic features which are also 
observed in the experiment for 12.5 Hz (see Fig. 6a in [10]). Beside the agree- 
ment in the generic features, there are three major quantitative differences. 
The period-2 state area between the harmonic response and the tongue is much 
thinner than in the experiment. An extended area of period A^ = 3 could not 
be found and the right limit of the period-2 state is too low compared to the 
experimental results. 

For a frequency of / = 0.2 (~ 16.4 Hz) the results are shown in Fig. 9. 
For a low susceptibility of x = 1-15 the peak starts to oscillate harmonically 
independently of the strength of AB. With increasing Bq the period-1 state is 
followed by a period of A^ = 4 for AB < 0.125. For AB > 0.125 the harmonic 
response is mainly replaced by the period-2 state, which is then replaced by 
the period 3. One notes the appearance of oscillations between two positive 
extrema inside the second regime for the low susceptibility case (Fig. 9a). For 
a high susceptibility of x = 2.5 the whole period diagram displays a band-like 
structure (Fig. 9b). For a fixed amplitude AB and increasing Bq, the period 
A^ = 1 appears first. Then either the period A^ = 6 follows for AB < 0.15 
or the periods N = 2 and A^ = 6 follow for AB > 0.15. The whole structure 
of periodic orbits ends with a broad band of period A^ = 5. This last novel 
feature is remarkable because no similar phenomenon has been observed in 
the experiment. For all tested frequencies in [10], the second regime gives way 
to the third regime by a period of A^ = 1 or A^ = 2. 

The comparison between the experimental and theoretical data generally shows 
a qualitative and partly a quantitative agreement with the dynamics of the 
peak. This agreement is achieved with a certain choice for the mass of the 
peak (17) and for the strength of the impact (19). The shown results are ro- 
bust against modifications of (17, 19) by a constant of 0(1). It is not necessary 
to fit parameters as the damping constant, the driving period, the critical field, 
and the resolution limit of the height in contrast to the minimal model in [10]. 
The other improvements are a more realistic nonlinear force term and the mul- 
tiplicative character of the driving. Our results at low and medium frequencies 
for X = 2.5 support the presumption that the MF used for the measurements 
of the dynamical behaviour has a higher susceptibility than given in [10]. The 
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experimental results at a high frequency of 23.5 Hz (see Fig. 7a in [10]) could 
not be found in our tested range of frequencies, 0.01 < / < 0.5. 



5 Summary 



In order to describe the complex and nonlinear dynamics of a single peak of 
the Rosensweig instability in an oscillatory magnetic field, we propose a model 
aiming at an analytical equation for the height of the peak at its centre. Our 
model approximates the peak by a half-ellipsoid atop a layer of magnetic fluid. 
By exploiting the Euler equation for magnetic fluids and the analytical results 
for a rotational ellipsoid, we obtain a non/mear equation for the dependence of 
the peak height on the applied induction (15). For static induction the quality 
of our proposed model is tested. It leads to the correct subcritical character 
of the bifurcation and gives the right width of the hysteresis compared with 
experimental results. 

For a time-dependent induction the effects of inertia and damping are in- 
corporated into equation (15). In correspondence with the experiments the 
dynamics is studied in a region, where the peak alternates between zero and a 
maximal height hmax- Our model shows not only qualitative agreement with 
the experimental results, as in the appearance of period doubling, trebling, 
and higher multiples of the driving period. Also a quantitative agreement is 
found for the parameter ranges of frequency and induction in which these 
phenomena occur. 

For low frequencies the response of the peak is harmonic for nearly any 
strength of the external excitation which is a superposition of a static part and 
an oscillatory part. The whole area of harmonic response is cone-like shaped 
in accordance with the experiment. For a medium frequency a structure of 
periods is found, where a tongue of high periodic and non-periodic oscillations 
appears. For low values of the amplitude of the alternating induction, the 
tongue directly follows the period- 1 state. For higher values of the amplitude 
the tongue is embedded in the period-2 state. The appearance and the loca- 
tion inside the parameter plane of an area of high periodic and non-periodic 
oscillations agree with the experimental data in the same frequency range. 

Beside the agreement with the generic features observed in the experiment 
at low and medium frequencies, the model predicts a novel phenomenon. For 
a frequency of about 16.4 Hz the peak oscillates with the period A^ = 5 
as the final period before the oscillations between zero and hmax end (see 
Fig. 9b). It would be challenging to seek a final period greater than 2 in the 
experiment, because for the studied frequencies in [10] the final oscillations 
have only periods of iV = 1 or A^ = 2. 
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In the dynamics of a magnetic fluid with a low susceptibihty a mixing of areas 
with different types of oscillations is found. For frequencies which are not too 
low, areas with oscillations between two positive extrema appear regularly 
inside areas with oscillations between zero and hmax- It would be interesting 
to test in experiments whether such a mixing can be observed for MF with 
low susceptibilities. 
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Fig. 1. Approximation of the MF peak by a half-elHpsoid with the vertical semiprin- 
cipal axis h and the horizontal semiprincipal axis R. The shape of the peak was 
measured in a setup used in [10] with the MF EMC 901 at B = 115.63 • 10"^ T 
(courtesy of A. Tiefenau). The height of the peak is ~ 8.2 mm. 
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Fig. 2. Sketch of an ellipsoid inscribed into the crest of a sinusoidal surface wave 
with the wave length 4 R. 
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Fig. 3. (a) Height of the peak h versus the strength of the static external induction i?o 
as solution of Eq. (15) for x = 0-01 (dashed line), x = 1-15 (long-dashed line), and 
2.5 (solid line). The critical induction for the subcritical bifurcation is independent 
of X) whereas the width of the hysteresis increases with increasing susceptibility, (b) 
The width of the hysteresis, Be — Bg, is plotted versus the susceptibility x of the 
magnetic fluid. 
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Fig. 4. The induction at the top of the peak Bh versus the height of of the peak h. 
Bfi {h) is scaled with respect to the external static induction Bq (the radius of the 
half-ellipsoid R). The solid line gives the solution of the approximation (10) with 
A = —1/2 and the filled squares indicate the results of the fully three-dimensional 
calculations [18] (courtesy of Matthies). The induction at the tip of the peak is 
approximated with an accuracy of 1.7%. 
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Fig. 5. Three different regimes for the temporal behaviour of the height of the 
peak at constant frequency and constant amplitude of the alternating part of the 
applied induction. For small values of the static induction the height is zero (left), 
for moderate strengths the height oscillates between zero and a maximum (middle) , 
and for large values the height alternates between a minimum and maximum both 
larger than zero (right). Only the second regime is analysed in detail. 
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Fig. 6. Periods of the peak oscillations in dependence of the static induction Bq and 
the amplitude of the alternating induction AB at a small frequency of / = 0.03 
(~ 2.5 Hz) for the susceptibilities x = 1-15 (a) and x = 2.5 (b). The peak oscillates 
harmonically almost everywhere in the Bq-AB plane. The area of harmonic response 
is cone-like shaped. The limit to the left is given by AB = Be — Bq for Bq < Be 
(solid line) and the limit to the right is given by AB = Bq — Bc for Bq > Be (dashed 
line). Slight deviations from these features appear for x = 1-15 (a). The colour code 
for the periods 1 to 30 is given in the legend. Periods above 30 and non-period 
behaviour are displayed in grey. 
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Fig. 7. Periods of the peak oscillations in dependence of the static induction Bq and 
the amplitude of the alternating induction AB at a medium frequency of / = 0.1 
(~ 8.2 Hz) for the susceptibilities x = 1-15 (a) and x = 2.5 (b). (a) Two disjoint 
parts appear, where for smaller values of Bq the peaks oscillates harmonically and 
for higher values of Bq the period-2 state emerges, (b) For high amplitudes AB 
harmonic oscillations and period doubling are only present. For smaller amplitudes 
of AB a tongue of high periodic (2 < A^ < 19) and non-period oscillations appears. 
The solid line indicates AB = B^ — Bq for Bq < Be and the dashed line marks 
AB = Bq — Be for Bq > Be- The colour code for the periods 1 to 30 is given in the 
legend. Periods above 30 and non-period behaviour are displayed in grey. 
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Fig. 8. Oscillations of the height of the peak with a period N = 5 (a) and A^ = 16 
(b) at a driving frequency of / = 0.1 (~ 8.2 Hz) for a MF with a susceptibility of 
X = 2.5. The other parameters are: (a) AB = 0.15 and Bq = 0.74, (b) AB = 0.225 
and Bq = 0.74. Note the different scales at the axes. 
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Fig. 9. Periods of the peak oscillations in dependence of the static induction Bq and 
the amplitude of the alternating induction AB at a frequency of / = 0.2 (~ 16.4 
Hz) for the susceptibilities x = 1-15 (a) and x = 2.5 (b). (a) The peak starts to 
oscillate harmonically independent of the strength of AB. With increasing Bq the 
period-1 state is followed by a period of A'^ = 4 (small AB) or N = 2 and 3 (large 
AB). (b) The whole period diagram displays a band- like structure formed by areas 
of period A^ = 1, 2, 5, and 6. The solid line indicates AB = Be — Bq for Bq < Be- 
The colour code for the periods 1 to 30 is given in the legend. Periods above 30 and 
non-period behaviour are displayed in grey. 
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